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Abstract 

We present a new family of exact solutions of the Einstein equations that may be 
interpreted as representing the propagation of a pair of solitons, in the background 
of a plane-wave collision spacetime. The family is constructed through the Khan- 
Penrose procedure, as an extension of a known metric from an interaction region. 
The metric in the interaction region is obtained as a diagonal solitonic perturba- 
tion of Rindler's spacetime, applying the Belinskii and Zakharov Inverse Scattering 
Method (ISM), with two real poles and one complex pole and its complex conju- 
gate. We use in this solution a non-standard renormalization procedure, obtaining 
solutions that contain two more parameters than the analogous solution that results 
applying the standard ISM. 

We analyze the asymptotic behaviour of this family of solutions in the limit where 
the determinant of the two by two Killing part of the metric vanishes. We find that 
there exists a curvature singularity in this limit, except when the free parameters 
contained in the solutions satisfy certain relation, in which the new parameters 
introduced by the non-standard renormalization procedure play an essential role. 
When this condition is satisfied, it is possible to find a transformation to a coordinate 
system where the metric is regular in the limit indicated above, and we show that 
the resulting collision spacetimes contain in that region a Killing Cauchy horizon 
instead of a curvature singularity, as in the general case. Finally, we analytically 
extend this subfamily through the horizon, and we find a curvature singularity in 
this extension, that may be considered as the result of the perturbation introduced 
in the interaction region by to the presence and propagation of the two complex 
poles. 
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1 Introduction 



The global structure of a gravitational plane-wave collision spacetime, and, in particular, 
the singularities produced after the collision because of the strong mutual focusing, is a 
subject of permanent interest in the area of exact solutions in general relativity. The first 
exact vacuum solution of this type was obtained by Khan and Penrose [p]]. It describes the 
head-on collision of two parallel polarized^] impulsive plane waves. The final result of the 
collision is a curvature singularity that can not be avoided by any observer (in this sense 
we can say that the singularity is spacelike). This type of singularity was also obtained 
by Nutku and Halil in , where they describe the collision of two impulsive plane waves 
with arbitrary relative polarization. 

Taking a different starting point, Szekeres 0, |J analyzed the collision of plane waves 
as an initial value formulation, obtaining an exact solution corresponding to the collision 
of two parallel polarized but otherwise arbitrary plane waves. A large set of solutions 
was obtained by Chandrasekhar, Ferrari and Xanthopoulos, through the use of relations 
between the mathematical theory of black holes and that of plane-wave collisions |5|, || 0, 
H ||. The crucial point to obtain these solutions is that, in both cases, Einstein's equations 
are reducible to the same Ernst equation. In all these solutions, a curvature singularity 
develops as the result of the gravitational interaction, except in § where a Killing Cauchy 
horizon is obtained, instead of the curvature singularity. This horizon is a 3-dimensional 
hypersurface where the spacetime plane-symmetry is broken, due to the fact that one 
of the two spacelike Killing vectors becomes null on the horizon. Chandrasekhar finds 
an analytic extension through this horizon, verifying that the null Killing vector becomes 
timelike through the horizon. He also finds a timelike curvature singularity in the extended 
region. This kind of behaviour, that is, the development of a curvature singularity, or the 
break-down of the plane symmetry by the formation of a Killing Cauchy horizon, is the 



general outcome of an arbitrary plane-wave collision. This result was proved by Tipler fTC 
and his arguments suggest that these singularities are peculiar features of plane waves, 
because singularities are also the consequence of a collision of self-gravitating plane waves 
of other fields with arbitrary small energy density. The stability of the Killing Cauchy 
horizon upon variations keeping the plane symmetry [[TT|] , and the asymptotic behaviour of 
a parallel-polarized plane- wave collision solutions near the singularity [jTj|], were analyzed 
by Yurtsever. In this last work Yurtsever concludes that the metric near the singularity 
is asymptotic to a Kasner solution. 

As already mentioned, the solution found by Szekeres was obtained through an initial 
value formulation. This procedure, however, presents the difficulty of having to solve 
this initial value formulation in general relativity which, even in the highly symmetric 
cases considered in plane-wave collision analysis, is a non-trivial task. This is probably 
the reason why most of the research that has provided exact solutions to the problem has 
been based on the so-called Khan- Penrose construction (KPC) m this construction the 
initial value formulation is side-stepped, and the procedure acts as a generating technique 
which, starting from a known solution to Einstein's equations, not necessarily related to 

lr This is equivalent to saying that the 2x2 part of the metric corresponding to the plane symmetry 
is diagonal. 
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the collision problem but satisfying certain conditions [13| , provides a full solution to the 
problem. This simplicity in the construction has the inconvenient, however, that one does 
not know in advance what are the initial conditions for which one has found a solution. 
Rather, the final outcome of the collision is given, and one has to interpret the initial 
conditions through an analysis of the solution. 

Essentially, in the KPC the spacetime is divided into four regions (labeled I, II, III, IV), 
separated by null hypersurfaces containing the shock fronts, in such a way that IV is flat, 
II and III represent the plane waves before the collision and I gives the interaction region 
after the collision. A previously known solution provides the interaction region I, and 
then it is suitable extended to the other regions [13| . As we shall see below, the metric 
may have, in general, discontinuous derivatives on the hypersurfaces where the extension 
is carried out, and this brings in the possibility of Dirac's 5-type singularities in the 
curvature tensor on these hypersurfaces. Depending on whether some energy conditions 
are satisfied, these singularities may be interpreted as corresponding to a surface layer of 
null dust (massless particles). Nevertheless, they should be absent if we want to restrict to 
vacuum solutions. The corresponding conditions were found in |]13||, and will be described 



below, when we verify them in the new solution we present here. 

The starting point in the construction of our plane-wave collision spacetime is an 
already known solution of Einstein's equations flUf , obtained with the Belinskii-Zakharov 
Inverse Scattering Method (ISM) ||15|| . This method provides a way of obtaining solutions 
of Einstein's equations that can be regarded as solitonic perturbations of a given seed 
metric. These solitonic perturbations display some of the features found in solitons of 
other non-linear system such as the Korteweg-de Vries (KdV) equation. For instance, 
at a given time the perturbations is concentrated in certain region and travels with a 
well defined speed. Specifically, the solitonic solutions are characterized by the number of 
real or complex poles introduced as a part of the generation method [|15|]. Complex poles 
(which appear always together with their complex conjugate to guarantee a real metric) 
are associated with the type of solitonic perturbation mentioned above, while real ones are 
associated with singularities in the resulting metric (the nature of these singularities has 
been analyzed in [|16|], [|17]]). An important distinction regarding the solitonic solutions we 
use here, is that they were constructed with a non-standard renormalization procedure, 
given in ||14|| , and, therefore they contain two extra real free parameters per soliton, as 
compared with the standard procedure. These extra free parameters play a crucial role 
in the determination of the global properties of the collision spacetimes. 

In this paper, we take for region I a family of solutions obtained in [14], which are soli- 
tonic perturbations of a seed metric given by Rindler spacetime, containing two real and a 
pair of complex conjugate poles. This choice was taken because each real pole provides a 
singularity that will be interpreted as a shock front in the plane-wave collision spacetime, 
needed by the KPC (see below), and each complex pole gives extra non-trivial gravita- 
tional structure that contributes to the interaction after the collision. This extra structure 
and the new free parameters due to the non-standard renormalization procedure are the 
main new features of this family of solutions, as compared with earlier work of Ferrari 
et.al. For simplicity we first describe the resulting plane-wave collision spacetime 

with only two complex poles (and of course, the two real poles needed for the KPC), and 
then we generalize this spacetime to an arbitrary number of complex poles. The choice in 
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the seed metric is inspired in the already indicated asymptotic results obtained by Yurt 
sever 



12|j , who proved that in certain cases the plane- wave collision spacetimes evolves 
near the singularity (in the interaction region I) to a Kasner metric, of which the Rindler 
spacetime is a particular case. The particular Rindler case is very interesting because of 
the fact that the curvature singularity that characterizes a Kasner metric is changed to 
a Killing Cauchy horizon for the Rindler spacetime. It is then naturally interesting to 
see if the solitonic perturbations of Rindler's spacetime that can be interpreted as the 
propagation of solitons on a plane-wave collision spacetimes, maintain this horizon or if, 
instead, in all cases they lead to the development of a curvature singularity. 

In Section 2 we present the solitonic family of metrics, solutions of Einstein's equations, 
and we analyze the possibility of performing the KPC We find that this construction is 
possible only for certain values of the free parameters of the family. We also confirm 
that outside this range, the resulting spacetime cannot be interpreted as a collision space- 
time. In Section 3 we compute the asymptotic expression for the collision metric near 
the singularity in the interaction region I. With this asymptotic metric we calculate the 
Kretschmann scalar and we find a relation between the parameters that define the family, 
such that this scalar does not diverge. This result suggests the possibility of the pres- 
ence of a Killing Cauchy horizon instead of the curvature singularity. This suggestion 
is confirmed in Section 4, by finding a coordinate transformation which, for the cases 
mentioned above, lets us extend analytically the collision spacetime through the horizon. 
We also prove that a new curvature singularity is present in this analytically extended 
region, related to the presence and propagation of complex poles in the original solution. 



2 The plane-wave collision spacetime 

The starting point in the construction of a collision spacetime through the KPC, is a 
solution of Einstein's equations admitting two commuting spacelike Killing vectors. It is 
well known that the metric for this type of spacetimes can be written in the form 

ds 2 = f(u, v) du dv + g a b(u, v) dx a dx b (1) 

where indices a, b can take values 1, 2 while x a denotes the Killing coordinates, and u, v 
are null coordinates. We define the region I of the collision spacetime, or interaction 
region, as the set of events with u > 0, v > and the metric given by (|l|). The set of 
events given by u > 0, v < is region II, while u < 0, v > defines region III, and finally 
u < 0, v < corresponds to region IV. The coefficients of the metric in regions II, III and 
IV are defined in terms of metric in region I as follows 

fW(u,v) = /«(«,()) 9%\u,v) = g2(u,0) 

f( in \u,v) = /«(0 >V ) g!* b n \u,v) = gV>(p,v) 

f^(u,v) = /«(o,0) g2P(u,v) = g%(0,0) 

where and g *^ are the metric coefficients in region K = I, II, III, IV. The metric 
thus obtained is well defined and continuous in a neighbourhood of u = 0, v = 0. It 
can be checked that if the metric in region I corresponds to a vacuum solution, then the 
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metric in regions II, III, and IV is also a vacuum solution. In fact, in IV the metric is 
flat while in II and III the metric coefficients depend on only one null coordinate. The 
KPC is therefore interpreted as the head-on collision of two gravitational plane waves, 
with region I representing the interaction region after the collision. It is clear that the 
metric may have, in general, discontinuous derivatives on the hypersurfaces u = and 
v = 0, and this implies the possibility of Dirac's 5-type singularities on the curvature 
tensor, which should be absent if we want vacuum solutions. It can be checked that a 



sufficient condition to avoid this singular behaviour of the curvature tensor is [13] 



lim a )V = lim a >u = (2) 

where a = \J\g a b\ (i-e. the square root of the determinant of the Killing part of the metric) 
and subindices in (|2]) indicate differentiation. 

We choose as the metric that represents the interaction region an already known 
solution of Einstein's equations [14|, obtained by applying the ISM [15] to a seed metric 



given by the Rindler spacetime. We choose this spacetime as the background for the 
solitonic perturbation for the following reasons. First, it is of interest to search for collision 
spacetimes that do not develop curvature singularities, and, since a collision spacetime 
evolves near the singularity in the interaction region asymptotically to a Kasner spacetime 
T2|, we are rather naturally lead to study the perturbations of a Rindler metric, a special 



case of a Kasner spacetime, which does not develop a curvature singularity. Second, 
we consider a solitonic perturbation that incorporates also complex poles, in order to 
analyze the effect of this new structure and its propagation, over that of the already 
known solutions [18|], where only real poles were included. 



The chosen family of solution was obtained applying a non-standard renormalization 
procedure [[bf], that adds two extra real parameters, 5 and 7, to those already present in 



the standard Belinskii-Zakharov ISM. The expression for the metric in region I is, 

ds 2 = fit, z)(dz 2 - dt 2 ) + Gn(t, z) dx 2 + G 22 (t, z) dy 2 

where the Killing coordinates are now x and y, and the functions f(t,z), Gn(t,z) and 
G 22 (t, z) are given by 

G n (t,z) = t 2 ^II^J (3) 

G ^(t,z) = (4) 




MM) = {{W k - z) + e ky J(W k - zf - t^j (k = l---n) (6) 

where e k ± 1, n = 4, W\ = z±, W 2 = z 2 + iou, W3 = z 2 — iu and = z 3 , with u, 
zi, z 2 , Z3, and C real constants. We further assume that z\ < z 2 < z 3 to ensure that 
the perturbation produced by the complex poles can be observed inside the interaction 
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region. Summarizing, assuming that we keep the values of u, z±, Z2, z 3 , and C fixed, for 
each choice of the we have a well defined two-parameter family of solutions, (given by 
5 and 7), for the following ranges of the coordinates (which represent a triangle in (z,t) 
coordinates) 

— 00 < x, y < 00 
zi < z < z 3 (7) 
Z\ — z < t , z — z 3 < t and t < 0. 

With the choice of coordinates in (H)-® we have that a = t (this is always possible in 
metrics of the form ([!])). The constants Wk are directly related with the poles in the ISM 



TE\ . The complex conjugate pair of poles related to W 2 and W 3 introduce a non-trivial 
gravitational structure on the incoming plane waves, but do not give rise to a divergent 
behaviour on the metric components. On the other hand, the real poles related to W\ and 
Wi, besides their contribution to the non-trivial gravitational behaviour of the incoming 
plane waves, introduce also a divergent behaviour in the function /. This function is 
singular for t = Z\ — z and t = z — z 3 because of the vanishing factors (/i 2 — t 2 ) 7 and 
— t 2 ) 7 in its denominator. 

The KPC can be performed if there exist two null hypersurfaces (e.g. u = and v — 0) 
such that @ holds, or, written in our coordinate system, lim M ^ + t >u = and lim,;^ + t >v = 
0. Since any coordinate transformation from (z, t) to a pair of null coordinates (u, v) must 
be of the form t = F\(u) +F 2 (v), z = Fi(u) — F 2 (v), then (||) implies that lim n ^ + F% tU = 
and lim^ + ^2,« = 0- But the jacobian of this transformation is 2Fi^ u F 2jV , and, therefore, 
(0) implies that the coordinate transformation must be singular for u = and v = 
0. Therefore, since the metric in the null coordinates must be well behaved on these 
hypersurfaces (a requirement of the KPC), then it must be singular in the chart (z,t), on 
the same hypersurfaces. 

Thus, a necessary condition to construct a collision spacetime from ([|) - (^) is to 
choose two Wk, say W\ and W4, real, so that the metric is singular on the null hyper- 
surfaces t = z\ — z and t = z — z%, where the wave fronts will be located. The existence 
of these singularities is not a sufficient condition for the construction of a collision space- 
time, because we must also explicitly find the singular coordinate transformation to null 
coordinates that ensures (0). We can determine this coordinate transformation from the 
behaviour of the function / near any of these two null hypersurfaces. For example, near 
t = Z\ — z we have 



(t + z- Zl yf 2 / 2 

with fo(t, z) finite and different from zero near the hypersurface. Near the other hyper- 
surface we have an analogous expression. It can now be checked that the following is the 
desired coordinate transformation 

u = (z-z 1 +t) 1 /° , u E (0, (^3 — ^i) 1 /' 7 ] 

v = {z 3 -z + t) 1 '* , v G (0, (z 3 — Zi) 1 ^] { ) 



with 



a =\ TT^ > (u' 7 + v°)-(z 3 -z 1 )<0 , 7 2 <2 (9) 

1 — 7^/2 
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because under this transformation the metric is changed to 



, 2 2 f(u,v) , , _ , . , 2 t(u,v) 2 , 2 

as = — a - — - — auav + & n (u,v) ax + — — ; -ay 

Gh(«,v) 

and the divergent behaviour coming from the factors (/ii 2 — t 2 )~"' 2 and (fi^ 2 — t 2 ) -1 ' 2 in (|3|) 
- (H) is cancelled by the jacobian of the transformation, and the metric is well behaved 
on these hypersurfaces. Finally, to carry out the KPC, it is necessary to check that 
holds. It is easy to see that 

t = ±[(u° + v°)-(z 3 -z 1 )] 



so we have that 



lim t u = lim - u^ 1 = 7 2 < 2 (10) 

and the same condition for 7 is found on v = 0. Summarizing, we propose as the interac- 
tion region of our plane-wave collision spacetime, a family of metrics obtained as solitonic 
perturbation of Rindler's spacetime, which necessarily have two real poles to carry out 
the KPC. The construction is performed by obtaining explicitly new null coordinates u, v 
such that the wave fronts are the hypersurfaces u = and v — 0, and the metric has 
continuous first derivatives across these wave fronts. This construction, however, is not 
possible for the whole family of solitonic metrics, but is restricted to the subfamily with 
7 2 < 2. 

The rest of the solitonic family, for 7 2 > 2, cannot be considered as a collision space- 
time, because the regions t — z\ — z and t = z — z^, which we would like to interpret as the 
wave fronts are, in fact, at past null infinity, as can be seen from the following argument. 
First, we perform the following coordinate transformation] 2 ]. 

u=(z-z l + t)- 1 ^' , ue ((2; 3 -2; 1 )- 1 /-',oo) 

u = (z 3 -z + t)-V«' , DG^-z^.oo) 1 } 



with 



a 



f- 1 < 1 1 \ / 



7 2 /2-l ' K u 



- 7 + - 7 )-(z 3 -z l )<0 , 7 2 >2 (12) 



Notice that with u, and v, defined as in (|TT|), the hypersurfaces t — Z\— z and t = z — z 3 



correspond, respectively, to u = 00 and v = 00. The null vectors (d/du) a and (d/dv) a 
are, therefore, past directed. 

Next, we consider a future directed null geodesic T(u), with tangent vector T a = 
— (d/du) a , that contains an arbitrary point in region (0), whose new coordinates are 
(uq,Vq). The coordinate u is not an affine parameter for this geodesic, because 

1 V a 1 — 1 uu 1 



2 This coordinate transformation is only valid for 7 2 > 2, but it is easy to check that the same ideas 
are applicable in the case of 7 2 = 2, with a coordinate transformation of the form u = ln(z — z\ + t), 
v = ln(z3 — z + t); arriving at the same conclusion as in the case 7 2 > 2 
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Figure 1: The ratio (^2) (4)/ (^2) (2) i n the interaction region, for the parameters e± = e 2 = 
1, £4 = —1, uj = 0.1, 5 = and 7 = 1. 

but, solving 




we may construct an affine parameter A such that T a = (d/d\) a . For a metric of the form 
([D, in coordinates given by (0), we have 

Y u = In l — - — I 

and this implies that 

X{u) = -Lu^w du - 

We know that the function f(u,v) appearing in the metric has the form f(u,v) = 
fo(u,v) +1 ) near the hypersurface t = Z\ — Z, with f finite and different from zero 
in the hypersurface. This asymptotic behaviour of /, together with the fact that when 
7 2 > 2 the hypersurface t = Z\ — z corresponds to u = 00, means that the affine parameter 
of the geodesic on this hypersurface is 

rod 

A(oo) = - / fo(U,v ) dU = -00. 

Juo 

Thus, the geodesic reaches the hypersurface with infinite affine parameter implying that 
this hypersurface is at past null infinity and, therefore, the solitonic spacetime cannot 
be considered as a part of a collision spacetime. On the other hand, for 7 2 < 2, the 
hypersurface that defines the boundary between regions I and III, corresponds to u = 0, 
and a similar argument shows that in this case a past directed null geodesic reaches this 
boundary with a finite affine parameter. 
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Figure 2: The ratio (^2) (4)/ (^2) (2) in the interaction region, for the parameters e± = e 2 = 
64 = 1, uj — 0.1, 5 = and 7 = 1. 



Finally, we are interested in comparing the geometric structure that results from the 
presence of the complex soliton poles, with that obtained in the case where these complex 
solitons are absent This comparison can be done through the study of the Weyl 

scalars. The complexity of the metric coefficients makes this analysis very difficult, so 
we concentrate our attention on the Weyl scalar ^2, which is different from zero only 
in the interaction region. This Weyl scalar has the properties that it involves only first 
derivatives of the metric functions, and that it is the only non-zero Weyl scalar which is 
invariant under simultaneous rescaling of the tetrad vectors l a and n a . So, it is plausible 
to think that if ^2 is singular, the metric will be singular. 

The comparison is performed by introducing the following null tetrad, defined in the 
interaction region, 



r 



n 



m 



m 




and computing the corresponding expressions for the Weyl scalar (\I/2)(4) of the four- 
soliton metric and (\I/2) (2) of the corresponding two real soliton metric, that is, a metric 
with the same parameters as the former, but without the inclusion of the two complex 
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solitons. In Figures 1 and 2, we plot the ratio (^2) (4)/ (^2) (2), for some particular choices 
of the free parameters. We can see that the relative behaviour of these two of metrics 
depends in a crucial way on the choice of these free parameters. In Figure 1 we can see 
that ( x I / 2)(4) becomes very different from (\E f 2)(2) in the region t — > 0, while in Figure 2 
we see that this difference disappears in the same limit. The propagating nature of the 
soliton perturbations is also clear in both figures. 



3 Asymptotic behaviour 

Once we have constructed the plane-wave collision spacetime from the solitonic subfamily 
of Einstein's equations solutions, the next natural step is to analyze the resulting space- 
time, in particular, in the interaction region. It is well known JTDI that as the result of the 
strong mutual focusing of the waves, we have in this region the development of either a 
curvature singularity, or a breaking of the plane symmetry, with the creation of a Killing 
Cauchy horizon. In the solution we present here, this singular behaviour appears in the 
limit t — > 0~ in the chart (t,z). The exact explicit expressions for the coefficients of the 
metric in that region are very complicated. On this account, and to simplify the analysis, 
we consider appropriate Taylor series expansions, instead of the exact expressions. We 
use these expansions to compute the Kretschmann scalar, looking for sufficient conditions 
to ensure a curvature singularity. As we shall see, for some choices of the free parameters 
5, 7, e k , (k = 1---4), this scalar diverges in the limit t — > indicating a curvature 
singularity, but we find a finite limit for other choices of these parameters. We can show 
that for the latter spacetimes the metric approximates a Rindler form, suggesting that 
an appropriate coordinate transformation may be used to prove that the singular limit 
t — > 0~ in the chart (t, z) is only a coordinate singularity. In fact, using this coordinate 
transformation we find that the hypersurface corresponding to t = is indeed a Killing 
Cauchy horizon, and that the metric may be extended through this horizon. This exten- 
sion will be carried out in the next Section, while in this Section we construct the required 
coordinate transformation. 

We first calculate the expansion for the metric coefficients in the limit t — > 0~ and 
z = const. This may be done for arbitrary number of solitons n, without the restriction 
to n = 4. We begin our study expanding the coefficients \i k near t ~ 

e k t 2 



MM) 



{W k -z)+e ky /(W k -z) 



2J(W k - z) 



2yJ(W k -z)* if e k = 1 

(t 2 /2)/J(W k - zf if e k = -1 



with k = 1 • ■ ■ n. Here we have to consider different possibilities. If W k is complex, we 
have, 

Hk(t,z) ~ 

while for W k real, we have only two possibilities, say, W\ = z\ and W4 = z 3 . If we recall 
that Z\ < z < z 3 we obtain 

j / (t 2 /2)/(^i -z) if e x = 1 
" I 2(zx - z) if ex = -1 



lii(t,z) 
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/i 4 (t, z) 



2(^3 - z) 
(t 2 /2)/(z 3 



if £4 
if €4 



Defining B- k {z) = +2J(W k 



(k=£ 1, 4), = 2(*i-*) and£ 4 (z) = 2(z 3 - 



we may summarize all these results as 

MM) 



B k (z) 
t 2 /B k (z) 



Let us assume that we have chosen the e k such that there are m functions /J, k (t, z) which 
behave as t 2 / B k (z) when t — > 0~, z = const, and n — m such that fi k (t, z) — > B k (z) in the 
same limit. Then, the first order of the asymptotic expression for the metric coefficients 
in this limit is the following 



9(t,z) 



t 



2(5+m7) 



n 



j=m+l 



BkAz) 



n? =1 B ki (z) 



f(t,z) ~ Q f 2(5+m 7 )(5+m 7 -l) 

where C% is a constant. Defining the function 

g(z) 



-iB kj (z) 



11% 



UT=iB kl (z) 
iB k Azr 



(2(6+m7)-l) 7 



and the constant 
follows 



miB ki (z) 

(5 + 7717), we may rewrite the asymptotic form of the metric as 



ds 2 ~ C x t 2 ^- l) g(z)^- 1 ^ {-dt 2 + dz 2 ) + cfe 2 + t- 2 ^" 1 ) ^(z)" 1 rfy 2 . (13) 

With this expression for the metric, we obtain the following asymptotic expression for the 
Kretschmann scalar, 

320 2 (0- l) 2 ' 



RabcdR 



abed 



-1) + 1) 

Cl t W-Dj( 2 )2(2fl)i4 ' 



We conclude that if the parameters 5 and 7 are such that <ft is different from zero or one, 
then the corresponding plane-wave collision spacetimes develop a curvature singularity for 
t — > 0~. On the other hand, if <fi equals zero or one, we may prove that we have a coordinate 
singularity in the chart (t,z) when t — > 0~, by finding a coordinate transformation such 
that the metric is well behaved in the same limit. To obtain the appropriate coordinate 
transformation, we notice, from (|T3"D, that the asymptotic form of the metric corresponds 
to a Rindler spacetime. For example, for = we have, 



ds 2 



-1 



-dt 2 + t 2 dy 2 ) + g(z)" 1 dz 2 + g(z) dx 



where the constant C\ was removed by a rescaling of the coordinates. The singularity for 
t = corresponds here to the vanishing of the factor t 2 in the coefficient of dy 2 . But, if 
we introduce new coordinates r and 77, defined by 



T 



t cosh(y) 



V 



-t smh(y) 



(14) 
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we have, 



dt 2 + t 2 dy 2 = -dr 2 + drf 



(15) 



and in this new coordinates the asymptotic expression for the metric takes the following 
form 



and the singularity for t — 0, which in this chart corresponds to r = rj, is removed. A 
similar transformation may be applied if <fi = 1, with the change y <-> x, to prove regularity 
for t — when 0=1. We can also check that one of the spatial Killing vectors becomes 
null in this hypersurface. In the chart (t, z, x, y) the Killing vectors are 



and then £ a becomes null on the hypersurface r = ±77. In fact, it can be checked that 
this hypersurface is a Killing Cauchy horizon. This is because this hypersurface has no 
border (its generators are the integral lines of the null and spatial Killing vectors £ a and ( a 
respectively, and the spatial vector (d/dz) a with z varying up to the "fold singularity"); 
and this hypersurface is C _1 . The only non-trivial step is to check the completeness of 
the null geodesic generator. But it can be checked that the vector £ a is geodesic on the 
horizon, and since it is a Killing vector, then is complete. 

4 Analytic Extension 

We proved in the previous Section that, for certain values of the parameters 5 and 7, 
the corresponding plane-wave collision spacetimes develop a Killing Cauchy horizon after 
the collision, but the metric is otherwise regular. This behaviour was proved by writing 
the metric in the interaction region in appropriate Rindler coordinates (t,t],z,x), such 
that the interaction region corresponds to |r| > with r < 0, and the horizon is on 
the hypersurface t — —\r}\. It is then possible to extend the collision spacetime through 
the horizon, by analytically extending the definition of the metric coefficients, that is, 
defining this coefficients for all possible values of r and rj, where they are regular. 

The global structure of the spacetime extended in this form may be separated in four 
regions. The first region, which we call C±, corresponds to the original interaction region 
of the collision spacetime, that is |r| > \rj\ and r < 0; the second, C 2 , and third, C 3 , 
correspond to \r\ < \rj\ and rj < or i] > respectively; and the fourth, C 4 , is delimited 
by M > M and r > 0. 

However, the explicit expressions for the metric coefficients in this Rindler coordinates 
are quite involved, and so, to study the existence of singularities in any of those regions, 
it is convenient to perform appropriate coordinates transformations in each region. For 



ds 2 ~ g{z)~ l {-dr 2 + dr] 2 ) + giz)" 1 dz 2 + g(z) dx 2 




while in the chart (r, 77, z, y) they have the following expression 
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example in C4, it is convenient to make the transformation t 2 = r 2 — t] 2 , tanh(x) = —rj/r, 
which is the inverse of (|14]) . We obtain here the same functional expression for the metric 
coefficients as in (|3])-([5|), but now t is positive. We therefore see that the analytic extension 
of the plane-wave collision spacetime through the Killing Cauchy horizon involves more 
than merely matching the (t, z, x, y) charts for t < and t > 0, because we also have to 
consider the regions C 2 and C 3 . 

The analysis of the extended spacetime in the regions C 2 and C3, is simplified by the 
following coordinate transformation: r = tsinh(y), rj = tcosh(y). It can be checked that 
the metric coefficients obtained are the same as those in but with the change 

t 2 = —t 2 . This observation is useful in the analysis of the behaviour of metric coefficients 
in these regions. It can be checked that the coefficient f(t,z) is singular in the limit 
z = z 2 and t — > u, due to the vanishing of the factors z)) 2 + (t) 2 ) with {k = 2, 3). 

Therefore, this singular behaviour is directly related to the presence of a complex pole in 
the ISM. To decide whether this singularity corresponds to a curvature or a coordinate 
singularity, we again compute the Kretschmann scalar in the limit z = z 2 and t — > u. 
For technical reasons, it is easier to perform this calculation in coordinates £, £ given 
by t 2 = u; 2 (£ 2 + 1)(1 — £ 2 ) and z = co>££, where for simplicity we choose z 2 = 0. We 
calculate the asymptotic expression for the Kretschmann scalar in the limit £ = and 
£ — > 0~, which it is equivalent to the limit z = and t — > uj~, and we obtain the following 
expression 

R abcd R abcd = ( ] C 4(72 ' 2) [1 + 0(0] 




with F = / (f = 0, C = 0) and / (£, C) is defined as /(£, C) = /„(£, C)/(£ 2 + C 2 ) 7 ', and it 
can be checked that /o(£, C) is the non singular part of /(£, () in the limit £ — > 0, ( — ► 0. 
We conclude that all these collision spacetimes, as defined in Section 2, have a curvature 
singularity in the limit z = z 2 and t — > because, as we proved in Section 2, all these 
collision spacetimes satisfy •y 2 < 2. 



5 Conclusions 

We have presented a new family of plane-wave collision spacetime, constructed from an 
interaction region, through the Khan-Penrose procedure. The interaction region corre- 
sponds to a diagonal solitonic perturbation of Rindler's spacetime, obtained applying 
the ISM, with two real poles and a pair of complex conjugate poles. We further use 
a renormalization procedure developed in |TJ| to obtain a family of solutions with two 
more parameters than the analogous solution obtained with the standard ISM. We find 
conditions on the free parameter that characterize the family of solutions such that it is 
possible to perform the Khan-Penrose construction. We conclude that only a restricted 
subfamily (given by j 2 < 2), can be considered as a plane- wave collision spacetime, be- 
cause only in this case the hypersurfaces corresponding to the gravitational wave fronts 
can be reached by a congruence of past null geodesic, starting somewhere in the presumed 
interaction region, with finite affine parameter. For the rest of the family, these hyper- 
surfaces are located at past null infinity, precluding their interpretation as plane-wave 
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collision spacetimes. 

We also study the perturbation introduced by the two complex poles to the already 



known plane- wave collision solution corresponding only to two real poles ||18|| . This per- 
turbation is shown in Figures 1 and 2, where we plot the Weyl scalar ^2 of the four-soliton 
metric divided by the same Weyl scalar of the two-soliton metric, for different choices of 
the free parameters. 

Next, we analyze the asymptotic behaviour of this family of solutions in the limit where 
the determinant of the two by two Killing part of the metric vanishes. This analysis was 
carried out in the general case of a metric with an arbitrary number n > of solitons. We 
find that the Kretschmann scalar diverges, and therefore we have a curvature singularity in 
this limit, except when the free parameters of the family of solutions satisfy some relation, 
where the new parameters introduced by the non-standard renormalization procedure are 
essential. This condition also depends on the behaviour of the pole trajectory functions 
fik(t,z), (k = l---n), in this limit. When this condition is satisfied, we find that the 
Kretschmann scalar is well behaved, and studying the asymptotic expression of the metric 
coefficients, that resembles that of Rindler's spacetime, we showed that it is possible to 
find a coordinate transformation where the metric is regular in the limit mentioned above. 
We proved that this subfamily of collision spacetimes develops a Killing Cauchy horizon 
instead of a curvature singularity. 

Finally, we analytically extend this subfamily through the horizon. The extended 
spacetime obtained contains two stationary regions with a curvature singularity. This 
curvature singularity, which is timelike in the sense that can be avoided by an observer, 
may be interpreted as resulting from the perturbation introduced in the interaction region 
by the presence of the two complex poles. 
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